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Abstract 

Numerical simulations of the NLSE (or GPE) are presented demonstrating emission of short 
pulses of the matter (light) density formed in the course of tunneling in wave-guided light 
and/or trapped BEC. The phenomenon is observed under various conditions, for nonlinearities 
of different signs, zero nonlinearity included. We study, both numerically and analytically, pulsa- 
tions of matter (light) remaining within the trap and use the results in order to induce emission 
of sequential pulses by properly narrowing the trap. This allows us to propose a mechanism for 
a realization of Atom Pulse Laser. 
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1. Introduction 

The research on dynamics of systems described by Non Linear Schrodinger equation 
(NLSE), such as weakly interacting Bose - Einstein condensates (BEC) treated in the mean 
field approximation where it is known as Gross-Pitaevskii Equation (GPE), and light prop- 
agating in nonlinear refractive media, has generated many interesting results, most of which 
have to do with 'solitonic' behaviors of these systems. The realization of BEC in 1995 has 
made it possible to trace bright and dark matter wave solitons experimentally 1 2 3 4J , and 
enhanced further vast theoretical work on the subject. The similarity between the two phys- 
ical objects, optical systems and matter waves, both described by the similar equations, 
has made it possible to study both systems in parallel and shed light on one system by 
performing experiments on the other. Furthermore, the Atom Optics field has emerged, and 
very interesting applications of the optics principle implemented for matter waves were sug- 
gested and realized 5 6 7 8 9J. Interesting mechanisms for realization of pulsed Atom laser 
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were also suggested. Refs.[TU] and then [TT] propose to use spatial variation of the scattering 
length so that the trapped condensate is accelerated towards a new minimum thus having 
a large enough energy to escape from the trap as a single soliton. Partial outcoupling is 
then achieved through an additional local maximum which causes splitting in the attrac- 
tive interatomic interaction region, thus giving rise to a multiple emission. The authors also 
note that the velocities of the ejected solitons are constant and almost identical. Another 
mechanism was proposed in ref. [TJ] where an elongated, quasi Id repulsive condensate is 
subject to simultaneous variations of the harmonic trap, from attractive to repulsive, and of 
the interatomic interaction, from repulsive to attractive, resulting in self coherent solitonic 
pulses. A special attention is given there to the stability criteria, mainly for avoidable col- 
lapse in the realistic multidimensional case. In refs. [13|14j a periodic soliton array generator 
was introduced in the model that contained two tunnel- coupled parallel cigar-shaped traps, 
in which one trap is a BEC reservoir and the other trap is a lasing cavity. 

Theoretical studies on dynamics of such systems have been carried out both numerically, 
perhaps the most common is the Split Step Fourier (SSF) method, and analytically, where 
some thorough and fruitful investigation were carried out |15I16|17] . Tunneling in one and 
two dimensional systems in the WKB approximation (with and without vortices) was ad- 
dressed in the recent papers [18119] . Behavior of BEC under the action of a time dependent 
field and tunneling through chaotic areas was analyzed in ref. [20]. The hydrodynamic for- 
mulation of the Schrodinger equation was originally proposed in ref. [21] and has been in 
use since both in quantum mechanics and optics (see, e.g. [22 23 24 25 26 27 ). An ana- 
lytical study of macroscopic tunneling dynamics within the hydrodynamic representation 
of NLSE has been carried in Ref. [28129] , This study employs the adiabatic approxima- 
tion making use of the two characteristic time scales typical of the macroscopic tunneling 
processes 30 31 32 33], and results in an almost complete solution for the dynamics of tunnel- 
ing of weakly interacting BEC. The most important feature of this solution is the observation 
of a tendency to formation of a dispersive shock wave in the outskirts of the potential trap, 
which propagates later on as a blip or pulse in the matter density or light intensity and which 
under certain conditions may form a bright soliton. It is important to emphasize that such 
pulses appear irrespective of the nonlinearity in the NLSE, and the show up, in particular, 
in the case of zero interaction (or zero Kerr nonlinearity in Optics). 

Approaching BEC as a quantum fluid has recently given shock wave dynamic phenom- 
ena its rightfully deserved attention. Dispersive shock waves and related behaviors in BEC 
dynamics have been predicted theoretically and observed experimentally. In classical com- 
pressible gas and fluid dynamics, shock waves are known to be traveling fronts of steep 
gradients. However, since GPE admits no dissipation effects, the concept of dispersive shock 
waves has been introduced. 34 35 36 37 38] These are waves that are generated and main- 
tained by dispersion instead of dissipation and are believed to consist of a number of spatial 
dark solitons and exhibit large amplitude oscillations. 

In this work we further investigate the dynamics of 1-d macroscopic tunneling by simu- 
lating temporary and spatial evolution of tunneling matter wave and light described within 
the frame work of NLSE/GPE. We carry out detailed simulations which demonstrate an 
agreement with formation and propagation of a blip [29\ regardless of the type of nonlin- 
earity. However, the latter causes differences in the further behavior of the outgoing pulse 
in the three cases, negative, zero and positive nonlinearity, which correspond to focusing, 
defocusing and 'nonfocusing' media in Optics, and respectively attractive, repulsive and zero 
interactions in BEC. 
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We also demonstrate how different initial states rapidly transform into the state optimal 
for the blip formation used in Ref. [29]. This indicates that the process of blip formation 
is not really sensitive to the initial conditions. We then concentrate on the matter (light) 
which remains trapped after the emergence of the blip, which is seen, both numerically and 
theoretically, to pulsate with the frequency mainly determined by the trap. 

In the last part of this work we merge the two studied effects, i.e. single blip emission 
and periodic pulsations, to search for conditions under which blips may be repeatedly and 
controllably emitted. Such 'blip train' might pave a road for realization of Atom Pulsed 
Laser. ( fl0|ll|12|13|14j ) We show how each parameter of such train is controllable. 



2. Simulation 



Our simulation program is based on the Interaction Pictures Runge-Kutta-4 method [39 
for solving the dimcnsionless GPE 
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where ^(a^i) is the normalized to one wave function of the condensate with N particles in 
it. The external potential 
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is a symmetric tunnel trap (see Fig.[T]) parameterized by the widths a and b. C is a constant 
which detemines the height Uq = max[U ex t{x)] ~ 1 fo the trap. Ui n t is the effective interac- 
tion potential which may be positive (repulsive) or negative (attractive), or even zero. Our 
calculations use the values N ^ int = ±0.15, 0. Then the initial normalized Gaussian packet, 
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whose peak is situated at the center of the well at the time t = 0, evolves dynamically 
according to the GPE equation fl}. The parameters of the trap are C = 0.1 (so that Uq is 
close to one), a = 5, b = 1, and the width of the initial Gaussian input is w — 7. Fig. Q] 
presents a sketch showing the trap and the initial packet. 



2.1. Application to Optics 

In parallel, another simulation program, the Free Beam Propagation Method (FBPM) 
|40j . simulating optical pulses propagation in nonlinear media by solving the NLSE 

2 VA +( — ) - n - &n(x)A + i( — J — + f — j ■n-n NL \A\ 2 A = (3) 

was used. Here A is the electric field amplitude in the paraxial approximation. 

In order to get a full correspondence between Eq. ([3]) and the GPE J!}, we considered 
the time independent process in the two dimensional x, z space. The light beam propagates 
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Fig. 1. The tunnel trap (dotted line) and initial packet (solid line) for the parameters C = 0.1, a = 5, and 
6=1, and w = 7. 

along the z axis so that the z coordinate plays the role of time in the GPE. x is the 1-d 
traversal direction and corresponds to the space coordinate in the GPE. The 'potential trap' 
is actually a waveguide due to the variation of refractive index along the x axis which can be 
fabricated by correspondingly varying the dopant concentration. The deviation An of the 
refraction index from its global value has the form 



where we chose C = 2.28 for the calculations. An initial Gaussian pulse centered at x = 0, 
z = inside the trap is set to propagate in such a medium with the Kerr nonlinearity. The 
relevant medium (e.g., for AlGaAs) parameters are refractive index n = 3.33 at x — > oo, the 
Kerr nonlinearity ri2 = 3 • 10 16 ^|j, the maximal deviation of refractive index An max = 
0.001, the waveguide width W = 2000/zm, and its length L = 35,000/im. The initial input 
is a Gaussian shaped pulse of the width w = 300/im and intensities in the range 700 to 2000 
watt. The wavelength is A = 3/xm. 

3. Results 

3.1. Emergence of the Blip 

The simulation results for zero, negative and positive nonlinearities are presented in Fig. 
[2l We clearly see how a blip is created and propagates with a constant velocity independently 
of the nature of the nonlinearity, i.e. the blip shows up also in the case of zero or repulsive 
interaction in BEC, or zero or defocusing Kerr nonlinearity in optics. It surely exists in the 
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Fig. 2. Blip emergence and propagation. The figures show continuous top view onto the time evolution of 
BEC systems. From left to right: non focusing (zero interactions), focusing (attractive interactions) and 
de-focusing (repulsive interactions) cases . In the focusing case, the blip transforms into a bright soliton. It 
is also clear that the remaining trapped packet exhibit pulsations, a feature that will be discussed below. 



negative nonlinearity case (see Fig. [2]- focusing case), i.e. attractive interaction or focusing, 
in which case the blip under proper conditions may transform into a bright soliton. [25]. 

Another universal feature of this phenomenon is the velocity of the blip propagation. As 
can be seen in Fig. [2] it moves with a constant velocity (along straight lines). Moreover, 
the velocity remains the same for different nonlinearities. Our analysis in Ref. [29) yields 
equation 



Vblip = Vshock = \ (5) 

V m 

for the velocity of the dispersive shock wave propagation. Here Uo is the height of the po- 
tential barrier in the BEC problem. In optics it is translated into the dimensionless quantity 



vuip = y ■ (6) 

We see an excellent agreement between the theory and simulations such as v t heory = 1.4 and 
Vsimuiation = 1-3 in the BEC case and Vtheory = 2.45e~ 2 = 0.33 and v simulation = 2.57e~ 2 = 
0.35 in optics regardless of the type of nonlinearity. 

As was outlined in our previous paper [53] the blip results from the tendency to a shock 
formation on the outskirts of the trap. Fig.[2]shows how the focusing (attractive) nonlinearity 
leads to a conversion of the pulse width, controlled by both nonlinearity and dispersion, 
remains unchanged at least on the timescale of our calculations. A special attention to 
formation of a bright soliton from the emitted blip and its stability was given in our previous 
paper [29 . In the two other cases, the propagating blip starts exhibiting oscillatory tails 
typical of the dispersive shock waves 34 35 36 37 . These oscillations are visibly stronger for 
the defocusing case. The shock waves are governed by the dispersion and by defocusing 
media (repulsive interaction in BEC). 



3.2. Pulsations inside the trap 



An interesting result is that the initial packet wider than the trap very rapidly decays 
into a narrow packet loosely related to the ground state inside the trap with high shoulders 
outside the trap. This early time evolution takes place for all packets whose initial width is 
wider than the trap. 

As described above the initial packet emits a pair of blips and forms a rather narrow 
packet inside the potential well, which is rather close to the ground state of the well, obtained 
when disregarding the tunneling. This narrow packet is a state with a long life time which 
does not essentially decay during our simulation time. However, we observe well pronounced 
pulsations of the narrow packet with roughly doubled trap eigen frequency. These pulsations 
are seen in Figs.[2]and are shown in the 3-d plots in Fig. ((3]). We also see that in the focusing 
case, the oscillations have higher frequency and the latter increases as the nonlinearity 
increases, while in the defocusing case (repulsive interaction in BEC) the pulsations are of 
lower frequencies. 

In order to analyze the pulsations within the trap we continue along the lines of the hydro- 
dynamic approach whose application to the tunneling problem was presented in our previous 
paper[29j. The NLSE/GPE for the complex function can be rewritten as two equations for 
two real functions: the continuity equation 

p t (x,t) + [p(x,t)v{x,t)] x =0. (7) 

for the particle density distribution p(x,t) — \^(x,t)\ 2 and the Euler-type equation 
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for the velocity field v(x,t) — —V x (p(x,t) with f(x,t) being the phase of ^(x,t). 
We assume the harmonic approximation 

V ext (x) = — - — (9) 

for the shape of potential inside the trap and for the time being neglect the nonlinearity. The 
density of narrow packet formed after the emission of blips is assumed to have a Gaussian 
shape 



Po(x) = j= exp <^ - — 

where uq does not necessarily coincide with the amplitude of zero point oscillations u ex t = 
yjh/mv in the potential ([9|) . It rather corresponds to another harmonic potential with the 
frequency u> = Ti/mu^. 

We look for the solution of eqs. ([7]) and © also in the Gaussian form 

" w = Wf" (?§>) (10) 



and 



v(x,t) = x^ln f(t) (11) 
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Fig. 3. Pulsations inside the trap after emission of a blip. From top left: zero, positive (repulsive/defocusing) 
and negative (attracting/focusing) nonlincarity. 
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where u(t) = ito/(t) and /(O) = 1. It means that the time dependence of the solution is 
parameterized by the single function f(t). The functions (p~0|) and (jlip solve eq. ([7]) and 
when substituted into (JS]) for A = yield equation 

f" = -^f+co 2 r 3 . (12) 

for the function f(t). This is the equation of motion for a 'particle' with the unit mass and 
with the coordinate / in the 'effective potential' 



The function f(t) oscillates in this potential well between its initial value /o = /(0) = 1 
(we assume also that /'(0) = 0) and the point /i = lu/v found from energy conservation 
considerations. These two values of the function f(t) correspond to the minimal and maximal 
widths of the pulsating Gaussian packet. 

To find the full solution for f(t) we note that eq. (fT2"| can be mapped onto the equation 
of motion for a radially symmetric 2d harmonic oscillator with the frequency v in which / 
is its radial coordinate. Then eq. (|12p is obtained for the motion with the conserved angular 
momentum oj. 

The solution satisfying the above initial conditions is 



to 2 

f(t) = d cos 2 vt + — sin 2 vt. (14) 

which describes the dynamics of the width of the trapped packet. 

It is interesting to indicate that there is a more general analytical solution 

*>-m*(m)»(m) (15) 



where 



b \W))= c,+iC > me «{^W)) (16) 

with the same parameterizing function f(t) as in (|10| . The parameters C\ and C2 should 
be found from the normalization of the density distribution (|15p and its initial shape. The 
error function in (|16|) is odd in x which makes the distribution (| 1 5[) asymmetric unless we 
assume that one of the parameters C\ or C2 zeros. If C2 = we return to the Gaussian 
distribution considered above. If, however, C2 = we have a symmetric distribution with a 
minimum at x = similar to what we see at the maximum amplitude of the pulsations. 

In the presence of nonlinearity (A ^ 0) we cannot get an exact solution, however we may 
consider eqs. (| 10[) and as an approximate trial functions and find the energy of the 
pulsating packet in units Tiv, 

e d ro P = \f 2 + W eff {f) (17) 

where the first term in the r.h.s. stands for the kinetic energy of the pulsations. The second 
term 
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Fig. 4. The effective potential l|18j l f° r various interaction parameters: NX/u = 0.1 — upper curve, 
NX/v 2 =0 — middle curve, NX/u 2 = —0.1 — lower curve, uj/v = 1.5. The horizontal line shows the 
level of oscillations between the points fg and f\ in the case of A = 0. 

represents the new effective potential (instead of (|13p ). which is just the internal energy of 
the packet at the given radius fuo obtained by averaging the BEC hamiltonian |T]) with 
the above trial function. It means that the oscillations of the packet are governed by the 
equation of motion 

r = -sf+<sr* ■ XN 



PV2k 

Fig. [4] presents the effective potentials for three values of the dimensionless nonlinearity 
NX/v 2 . We see in this figure that the repulsive interaction decreases the minimal width and 
increases the maximal width, i.e. it increases the amplitude of pulsations. A more detailed 
analysis shows that the period of pulsations also increases — the frequency lowers. The 
attractive interaction acts in the opposite direction and the frequency grows. This type of 
behavior agrees with the simulation results which can be clearly seen from Fig. [3] - the 
number of peaks in the same time interval varies for the three graphs corresponding to the 
different values of the interaction parameter. 

The graph in fig. [3] for the nonfocusing case (zero interaction) shows that the period is 
roughly T = 8. To estimate the trap frequency, which is nearly harmonic near the center of 
the well, we fit it to a parabolic function from which we deduce the frequency to be around 
Vtrap — 0.0737. Therefore the expected period would be T = ~Kjv trav — 8-18 in a very good 
agreement with the theoretical prediction. 

The simulation results for the NLSE, eq. ©, show that the fraction of the initial pulse 
which stays trapped along the propagating axis, exhibits pulsations of period ~ 5000m fi. 
To discuss the frequency equivalent for light propagation we look at the second term in eq. 
(|3|). Taking only the inside part of the trap and removing the constant offset, this term can 
be rewritten as 
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~-(x) a - B - An — ilkY (19) 

and is to be fitted to the harmonic potential v 2 x 2 /2 for trapped particles. Then we get 
T = tt/v — 7200/im for 'the spatial period', which is close to the simulation result T S i m ~ 
5000/im. 

It is worth mentioning here that various type of oscillations have been discussed in lit- 
erature (see, review [16j ) as a rule using the Thomas- Fermi approximation which neglects 
the quantum pressure term. This approximation works especially well for small amplitude 
oscillations with the repulsive interaction when the density distribution is rather close to 
uniform. It is emphasized here that the density distribution in the tunneling problem we 
discuss here is strongly nonuniform and the pulsations discussed here are strong so that the 
role of the quantum pressure is of an utter importance. 

4. Atomic Pulsed Laser 

Combining the two discussed effects, i.e. blip emission due to macroscopic tunneling and 
in-trap periodic pulsations, one can propose a new mechanism of controllable multiple blip 
emission. Such effect may lead to a realization of Atomic Pulsed Laser. The mechanism we 
propose is based on narrowing the well to a desired width, according to eq. (| 13[) at times 
that match the pulsating packet being at its maximum width and allows for outcoupling 
matter waves pulses. In addition to the simple realization, a full control and modulation 
possibility of several solitonic features, also in the course of emission, is at hand. In order to 
demonstrate this we consider the potential 

(20) 

with the varying in time parameters 

)oi, C\ t <t\ 
a 2 ,C 2 t<t 2 (21) 



Variations of ait) and C(t) are fitted in such a way as to keep the height of the barrier 
unaltered if identical velocities of the pulses are desired. Alternatively they may be chosen 
to vary the barrier height in order to reach the desired variation of the velocities of outcoupled 
pulses. a(t) is the variable width of the well. The times t\, t 2 ... are set to match the period 
of pulsations. It is important to stress here that the latter changes according to the width 
of the narrowing trap as was proved earlier. Therefore, on-switch times alter accordingly. 

Three pairs of blips emitted from the trap are seen in fig.0 The first pair of blips is emitted 
at the start of the process and the two other pairs are emitted due the two narrowings at 
t = 8 and t = 13 with the widths a — 5, 4 and 0.7 respectively. Since the nonlinearity plays 
a minor role in the emission process this simulation is carried out without it. 

Velocity modulation. As stated above, our analysis, based on the hydrodynamic ap- 
proach, predicts a one to one relation between the velocity and the height of the barriers, 
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Fig. 5. Emission of the second and third pairs of blips by the trap twice subsequently squeezed at the maximal 
width of pulsations. 

leading to a viable control of the outcoupled pulse velocities also in the course of sequential 
emissions. In this case the velocity of the blips is altered according to 



Fig. [6] shows the velocity variation for a single pair of pulses during its propagation, and for 
two pairs of pulses which are made to propagate with different velocities. 

The velocity difference between the two consecutive outcoupled pulses has been created by 
temporally changing the barrier height. The restriction here being that the first propagating 
outcoupled pulse has to be far enough from the barrier region, so that it will not be effected 
as well. 

Temporal interval modulation and control Temporal spacing modulation can be 
achieved by considering different numbers of pulsation periods between two emissions. An- 
other feature that should be taken into consideration is the time of the blip formation. This 
time is shown in ref. [29j to vary according to the trap parameters as 




vuips = \/2C(t)s(i) 
where s(t) is the solution of the transcendental equation 
tanh (a(t) s) 



(22) 



(23) 



a(t) 





(24) 



Therefore the formation time becomes shorter for the narrower barriers. 
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t 

Fig. 6. Velocity modulation achieved by raising the barrier height. Left: modulation of single pulse velocity, 
during its propagation. At t = 15 the barrier peak is raised from ~ 1 to ~ 4 Right: production of two 
consecutive pairs of pulses with different velocities. At t = 15 the trap is narrowed from a=5 to a=2.5 and 
the barrier peak is raised from ~ 1 to ~ 2 

Stability The pulse creation is not related to the nonlinearity of the GPE and the present 
discussion without loss of generality is limited to the zero interaction. However, it is of course 
evident that for an operating pulsed laser the pulses have to transform into bright solitons 
and moreover need to be kept stable and not disperse. This can be achieved by turning 
on attractive interaction along the way with lower limit according to the soliton formation 
condition 

U mt • |V>(0, T)\ 2 ~ 6.5a- 2 + 0.03C/ (25) 

found in ref. [29] for the shape of trap potential which is rather close to the trap chosen 
here. In fact the parameter values in eq. (|25p may slightly change with the shape of the 
particular trap. Here T is the on-set time. 1-0(0, T)\ 2 the density in the center of the trap at 
t = T. If the condition [25] is satisfied, the emerging pulse may transform into a stable bright 
soliton. In 2d and 3d case some stricter conditions as say Vakhitov-Kolokolov 41] should be 
satisfied. Since the emerging pulses are much more dilute than the initial packet, stability 
against collapse in 2d and 3d case is easier achievable. The masses and widths of the pulses 
can be also found as described for the single emission event in ref. [29"]. 

In principle we can consider restoring the original width of the trap after each emission 
and repeat this process as many times as necessary. This mechanism will work as long as 
there is matter left within the trap and its coherency is held. 

5. Conclusions 

Certain aspects of the universality of the emission of blips in the course of macroscopic 
tunneling were asserted, using two simulation programs, one of which applying directly to 
Optics. Pulsations of the trapped packet after emergence of the blips were shown by the 
hydrodynamic approach and by numerics to occur with the period matching twice the trap 
frequency, which inspired us to suggest a mechanism for realization of Atom Pulsed Laser. 
The blips of the matter density may be emitted by controllably narrowing the trap. The 
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velocity of the outgoing pulses may be controlled by varying the height of the barrier. We 
show also how to control the times of the pulse emissions and their masses. 
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